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ABSTRACT 

We present Doppler imaging and Balmer line analysis of the weak-line T Tauri star 
TWA 6. Using this data we have made one of the first attempts to measure differential 
rotation in a T Tauri star, and the first detection of a slingshot prominence in such 
a star. We also show the most direct evidence to date of the existence of solar- type 
plages in a star other than the Sun. 

Observations were made over six nights: ll-13th February 2006 and 18-20th 
February 2006, when spectra were taken with the UCL Echelle Spectrograph on the 
3.9-m Anglo-Australian Telescope. Using least-squares deconvolution to improve the 
effective signal-to-noise ratio we produced two Doppler maps. These show similar fea- 
tures to maps of other rapidly rotating T Tauri stars, i.e. a polar spot with more spots 
extending out of it down to equator. Comparison of the two maps was carried out 
to measure the differential rotation. Cross-correlation and parameter fitting indicates 
that TWA 6 does not have detectable differential rotation. 

The Balmer emission of the star was studied. The mean Ha profile has a narrow 
component consistent with rotational broadening and a broad component extending 
out to ±250 km s~^. The variability in Ha suggests that the chromosphere has active 
regions that are cospatial with the spots in the photosphere, similar to the 'plages' 
observed on the Sun. In addition the star has at least one slingshot prominence 3i?* 
above the surface - the first such detection in a T Tauri star. 

Key words: stars: pre-main-sequence - stars: late-type - stars: rotation - stars: spots 
- stars: individual: TWA 6 - stars: magnetic fields 



1 INTRODUCTION 
1.1 T Tauri Stars 

T Tauri stars are pre-main sequence equivalents of Sun-like 
stars. They are youn g, have low ma sses (0.3 - 1.0 Mq) and 
spectral types F-M (|Havashi|[l966l ). The T Tauri phase is 
considered to begin when the star becomes optically visible 
and continues until hydrogen burning begins and the star 
has reached the post-T Tauri phase at an age of the order 
of 10^ years (depending on spectral type). 

Very young stars fall on Hayashi tracks on the 
Hertzsprung-Russell diagram, i.e. vertical tracks of constant 
temperature. As the star approaches the zero-age main se- 
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quence (ZAMS) it can either continue to follow the Hayashi 
track all the way onto the ZAMS (and remain fully convec- 
tive), or develop a radiative core and turn onto the Henyey 
track and evolve onto the ZAMS with constant luminos- 
ity and increasing temperature. Which of these evolutionary 
paths the T Tauri star follows is thought to depend on the 
mass of the star but observations of TTS with ages close to 
this turn-off will allow us to characterise this stage of stellar 
evolution more fully. 

T Tauri s tars fall into two subcategories - classical 
and weak- line j BertoutI 1 19891 and references therein) . Ob- 
servations of classical TTS (CTTS) show excess ultravio- 
let and infrared emission and strong emission line activity, 
characteristics which are due to the presence of a circum- 
stellar disc. The star-disc interactions lead to strong, ir- 
regular variability that is predominately due to accretion 
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of material from the disc onto the surface. The disc has a 
braking effect on the star, causing C TTS to rotate more 
slowly than their discle ss equivalents (jBouvier et al ] |l993l : 
ICollier Cameron fc Campbell. 19931 ). 

Weak-line T Tauri stars (WTTS) have weaker emission- 
line activity and show less evidence of accre tion. WTTS ap- 
pear to be older than CTTS on average (jArmitage et al.l 
I2OO3I ). suggesting that when TTS form they have circum- 
stellar discs which dissipate as the star evolve s, but disc life- 
times are variable, ranging from 0.1-10 Myr (|Bouvier et al.l 
liggTl ). WTTS are variable, but in contrast to CTTS the vari- 
ation is usually periodic as it is due to cool spots that go in 
and out of view as the star rotates. These cool spots form 
when magnetic flux tub es emerge through the photosphere. 

Doppler imaging (|Vogt fc PenrodI Il983l ) maps spots 
on stellar photospheres. Doppler imaging of rapidly rotat- 
ing, magnetically active stars, such as WTTS, often finds 
spots at high latitudes a nd even at t he po les, see e.g. 
I Rice fc Strassmeieij l|l996l ). IStrassmeieij (|2002l ). in contrast 
to the Sun where they are normally within ±30° of the equa- 
tor. One of the favoured explanations is that Coriolis forces 
in rapidly rotati ng stars cause the flux tubes to emerge at 
highe r latitudes jSchuessler fc Solankil Il992l : ICaligari et al.l 
1 19941 ). 

In the past Doppler imaging of cool stars has been 
limited by its challenging signal-to-noise requirem ents. The 
techn ique of least-squares deconvolution (LSD. iDonati et al.l 
1 19971 ) has made it possible to produce Doppler images of 
stars for which the signal-to-noise ratio (SNR) would be too 
low if individual spectral lines were used. In LSD the under- 
lying rotationally broadened line proflle for each spectrum 
is found by deconvolving a line list from the spectrum. This 
makes it possible for late-K T Tauri stars to be imaged, 
where previously the SNR requirements would have been 
prohibitive. 

Sunspots are accompanied by prominences, i.e. loops 
of material that follow the magnetic fleld lines emerging 
through the surface of the star. High-lying prominences, 
known as 'slingshot prominences', typically appear as fast- 
moving absorption transients and have bee n observed in a 
number of young stars such as Speedy Mic (iDunstone et ahl 
20061) and AB D or (jCollier Cameron fc RobinsonI 1 19891 : 



Donati et al] 1 19991 ). but not in any TTS to date. These 



prominences are viewed as a tracer for large-scale magnetic 
fields. If a T Tauri star was shown to be able to support 
such high-lying prominences for several rotational periods it 
would suggest the presence of a large-scale field that does 
not decay away quickly outside the stellar radius. 



1.2 Differential Rotation in Stars with Deep 
Convection Zones 

Fully convective stars are believed to have a distributed mag- 
netic dynamo {a^), arising from turbulent motions or Corio- 
lis forces, rather than a solar-type dynamo {aQ.) that is due 
to both convection and differential rotation, as seen in stars 
with radiative cores. 

Differential rotation, when rotational velocity is a func- 
tion of latitude, is expecte d to be less apparent with incr eas- 
ing convection zone depth (|Kitchatinov fc Rudigei|[r999l ) , so 
stars with larger convection zones are expected to rotate 
more like a solid body. Differential rotation in the Sun gives 



rise to the S7 effect that produces a large-scale toroidal field 
from the small-scale poloidal field. Most models for an o? dy- 
namo rule out the possibility of there being differential rota- 
tion in f ully convective stars (e.g. iChabrier fc Kiiker 2006i ). 
However iDobler et al.l (|2006l ) suggest fully convective stars 
have significant differential rotation, although this would not 
be Sun-like, and takes the form of a meridional circulation. 

Observations to date have only attempted to detect 
the latitudinal differential rotation seen in the Sun. A re- 
cent measurement for t he fuUy-convective M4 dwarf V374 
Peg (|Morin et al.1 12007| ) found a weak differential rotation 
(around 0.1 of the solar value). Differential rotation has been 
observed to fall with d ecreasing phot osphcric temperature. 
For example, figure 2 in lBarnes et al. I {2005a,) plots differen- 
tial rotation against temperature. The trend seems to imply 
that differential rotation approaches zero below 4000K but it 
is uncertain where it becomes undetectable as there has been 
no measurement to date of stars with temperatures between 
3800-4500K. In addition there have not been any reliable 
measurements of differential rotation in T Tauri stars. Such 
measurements are needed to confirm that TTS follow these 
trends too. 



1.3 TW Hydrae Association 

Nearby young co-moving groups are a rich source of T Tauri 
stars for observation. Our chosen targets lie in the TW Hy- 
drae association (TWA) , a group of pre-main sequence stars 
lying at a distanc e of abou t 51 ± 5pc and with an age of 
8 Myr ({Zuckerman fc Sond [20041 ). Of the 24 likely mem- 
bers of the TWA , 16 ha ve had their periods measured in 
iLawson fc Cr^ |200i). The age of this association is an 
important one as it is corresponds to the mean age of cir- 
cumstellar discs and hence the time when T Tauri stars can 
begin to spin up, so the members of this association are inter- 
esting objects for study. The star that gives the association 
its name (TW Hya) is a classical T Tauri star, spectral type 
K7. Its disc has been directly imaged i nlKrist et ahl (I2OOOI) 
and its magnetic field was detected in I Yang et al.l (|2007h . 
Although it is rotating fairly rapidly (P = 2.8d) it is almost 
pole-on and is therefore not a suitable candidate for Doppler 
imaging. 

This paper concerns another member of the TW Hy- 
drae association - TWA 6. It has a spectral type of K7 
(jWebb et al.ll 19991 ). correspo nding to a temperatur e of about 
40 00K for a T Tau ri star l|Cohen fc Kuhil Il979l ). Models 
bv lSiess et all (|2000l ) indicate that it is fully convective or 
just beginning to develop a radiative core. Photometry has 
revealed a large lightcurve am plitude, with a variatio n of 
0.49 magnitudes in the V band llLawspn fc Crause |2005| ). Its 
rapid rotation, v sin i = 55 km s~^ (IWebb et al.lll999l ) and 
period of 0.54 days (|Lawson fc Crausell2005l ) make it an ex- 



cellent target for Doppler imaging. In this paper we present 
Doppler images of TWA 6, and a discussion of whether it has 
differential rotation. We also analyse the Balmer line emis- 
sion in order to study the behaviour of the chromosphere. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were carried out at the Anglo- Australian Tele- 
scope (AAT) using the University College London Echelle 
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Spectrograph (UCLES) and the EEV2 2k x 4k detector be- 
tween ll-20th February 2006. In order to detect differential 
rotation observations were made with a gap of 4 nights. By 
doing this two maps can be made of the star allowing differ- 
ential rotation to be det ected, or an upper limit placed on its 
value ijPetit et al.ll2004l ). With three nights of observations 
on either side of this gap full phase coverage was achievable 
on these two separate occasions. 

We took spectra of the targets, telluric and spectral 
standard stars, and calibration frames (see tables [l] and [2]) . 
UCLES was centred at 5500A giving a wavelength coverage 
of 4400 - 7200 A. A slit width of 1.0" was used giving a 
spectral resolution of 45000. 

The weather on the first three nights was reasonably 
good, and although conditions and seeing were variable (see- 
ing ranged from 1.2-2.2") in total only 3 hrs were lost. 36 
good quality spectra of TWA 6 were obtained, as well as 
27 spectra of our secondary target TWA 17 (Skelly et al, in 
prep.). On the second set of nights the weather was of vari- 
able quality and more time was lost. Of a potential 24 hours 
we had 13 hours of observing time, with seeing between 1.3 
- 2.2". Most of this time were concentrated on TWA 6 and 
a further 39 spectra were obtained. An additional 1.5hr of 
observing time was obtained on 17th Feb 2006 due to the 
early finish of aluminising the telescope mirror. 
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Figure 1. The dashed line shows the mean deconvolved profile 
for the first dataset (11— 13th Feb 2006). There is an absorption 
feature due to sunlight reflected off the Moon at km s~^. The 
solid line is the mean of the profiles after this feature was removed 
and a slope was removed. Mean errors are shown for the data after 
the feature was removed. The errors close to km s~^have been 
increased due to the subtraction. 



2.1 Data Reduction 

Data reduction was c arried out wit h the STARLINK rou- 
tine ECHOMOP (Mills et al.lll997l '). Frames were bias sub- 
tracted and flatfielded, cosmic rays were removed and the 
spectra were extracted. Sky subtraction was carried out by 
subtracting a polynomial fitted to the inter-order gaps. The 
spectra were normalised using a continuum fit to a star of 
the same spectral type. The spectral type (K7) was verified 
by eye by identifying the best temperature fit to the spec- 
tra in comparison with the template spectra. Wavelength 
calibration was carried out using ThAr arc spectra. 



3 LEAST-SQUARES DECONVOLUTION 

As the 36+39 spectra did not have a high enough SNR for 
Doppler imaging we used least-squares deconvolution (LSD, 
described in Donati et al...l997i ) to extract the underlying ro- 
tationally broadened profiles. LSD uses matrix inversion to 
deconvolve a line list from a spectrum to leave a high signal- 
to-noise pro file. The program SPDECON (implemented as 
described in I Barnes" et al.lll998l ) was used to do this. We used 
a line list from th e Vienna Atomic Line Database (VALD, 
iKupka et aP l2000[ l for a 4000K star. Over 4000 images of 
2000 spectral lines were used in the deconvolution leading 
to a signal gain of ~ 60. For the deconvolution, the lines 
were weighted by counts^ /variance. The effect of this is to 
favour regions of the spectrum with higher SNR, and the 
centres of the orders rather than the wings. The resultant 
line profiles had a mean centroidal wavelength of 5530A. 

There was a small consistent tilt in the profiles which 
was removed by fitting a straight line between the two sec- 
tions of continuum on the left and right sides. This tilt 
changed depending on which spectral standard was used and 
hence is most likely due to a small mismatch between the 



temperatures of the star and the spectral standards. In ad- 
dition there was a feature near km s"'^ in many of the 
profiles. It was identified as being due to moonlight in the 
spectra as it only appeared to be present when the Moon 
was above the horizon (Fig[l}. 

As we did not have a solar spectrum we deconvolved the 
hottest template spectrum (K3) we had taken with UCLES 
to model the shape of this feature. We were unable to find 
a pattern to its depth as it is difficult to disentangle from 
variations due to surface features. We failed to find a re- 
lationship between its equivalent width and the zenith dis- 
tance, time or seeing, although the absorption signal was 
stronger during bad weather. Hence for all spectra taken 
when the Moon was up, a constant shape was subtracted 
from all profiles taken within the same night (see Fig. [TJ. 
This was calculated by inspecting the mean of all profiles, 
and subtracting a profile such that the absorption feature 
was not obvious in the mean. A scale factor was calculated 
separately for each night. The errors for the pixels affected 
were increased so that these pixels would have a lower influ- 
ence on the final image produced. This was done by adding 
in quadrature the error produced by SPDECON and the 
amount that was subtracted from each pixel. 

Each of the deconvolved profiles has been normalised 
by dividing through by the mean profile. In Fig. [2] the nor- 
malised profiles have been stacked and shown in greyscale. 
As spots appear in emission the lighter regions in these im- 
ages show the spots being carried across the stellar disc. Spot 
groups are apparent at the same phases in both, particularly 
around phases 0.05 and 0.55, corresponding to longitudes of 
340° and 160°. This suggests that the spot distribution is 
similar at both epochs. 
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Table 1. Table setting out the observations made on the nights of 11/02/06 to 13/02/06. 
Spectral types in the 'Comments' column refer to spectral standards. The signal-to-noise 
ratio is the maximum for the set of exposures. 



V_7 U J ci^L 


Date 


TIT Start 


Exp. 

time [s] 


No. of 

Frames 


S /N 


VJtJlillllCllliB 


v.yUiiu.i Liuiia 


HR 3037 


11 


09:42:56 


20 


1 


98 


Telluric 




HR 3037 


11 


09:46:15 


120 


1 


164 


Telluric 




HD 34673 


11 


09:58:05 


200 


1 


68 


K3 




TWA 6 


11 


10:25:43 


900 


6 


53 


Target 1 




HR 3037 


11 


12:10:19 


20 


1 


105 


Telluric 




HD 52919 


11 


12:16:17 


200 


1 


79 


K5 




TWA 6 


11 


12:22:44 


900 


5 


63 


Target 1 




SAO 116640 


11 


13:45:02 


450 


1 


110 


K7 




TWA 6 


11 


13:55:31 


900 


5 


58 


Target 1 




TWA 6 


11 


15:17:37 


900 


5 


62 


Target 1 




TWA 17 


11 


16:43:21 


1200 


1 


48 


Target 2 




TWA 6 






900 


4 


53 






GJ 204 


12 


09:40:18 


120 


2 


58 


K5 




GJ 3331A 


12 


10:01:49 


900 


1 


16 


M2 


Diffuse cloud 


TWA 6 


12 


10:23:33 


900 


1 


19 


Target 1 


Diffuse cloud 


TWA 6 


12 


10:42:47 


900 


1 


3 


Target 1 


Aborted 148s 


HD 19007 


12 


11:04:33 


200 


1 


40 


K4 


Heavy cloud 


TWA 17 


12 


12:42:06 


1200 


1 


31 


Target 2 




HR 3037 


12 


13:06:28 


20 


1 


90 


Telluric 




TWA 17 


12 


13:13:37 


1200 


2 


33 


Target 2 




TWA 17 


12 


14:06:49 


1200 


3 


42 


Target 2 




HR 3037 


12 


15:26:56 


20 


1 


40 


Telluric 


Cloud 


TWA 17 


12 


15:42:30 


1200 


2 


28 


Target 2 


Cloud 


HR 3037 


12 


16:48:51 


20 


1 


87 


Telluric 




TWA 17 


12 


16:55:16 


1200 


3 


36 


Target 2 


Some fog 


HR 3037 


13 


09:32:42 


20 


1 


49 


Telluric 




HD 52919 


13 


09:35:50 


200 


1 


90 


K5 




SAO 116640 


13 


09:40:57 


400 


1 


82 


K7 




GJ 3331A 


13 


09:50:39 


900 


1 


18 


M2 




TWA 6 


13 


10:09:24 


900 


9 


37 


Target 1 


Seeing > 2" 


GJ 204 


13 


12:36:26 


120 


1 


67 


K5 




TWA 17 


13 


12:44:30 


1200 


2 


23 


Target 2 


Seeing > 2" 


TWA 6 


13 


13:31:51 


900 


1 


39 


Target 1 


Seeing > 2" 


TWA 17 


13 


13:49:34 


1200 


2 


30 


Target 2 


Seeing > 2" 


TWA 6 


13 


14:33:20 


900 


1 


42 


Target 1 


Seeing > 2" 


TWA 17 


13 


15:01:47 


1200 


2 


27 


Target 2 




TWA 6 


13 


15:35:40 


900 


1 


44 


Target 1 




TWA 17 


13 


15:54:08 


1200 


4 


31 


Target 2 




TWA 6 


13 


17:22:43 


900 


1 


49 


Target 1 




TWA 17 


13 


17:40:59 


1200 


4 


34 


Target 2 





4 DOPPLER IMAGING 
4.1 Doppler Imaging Code 

The unnormalised profiles and corresponding errors pro- 
duced by SPDECON are us ed as inputs for the program 
DoTS (IColher Cameronll200ll ). The program calculates the 
surface temperature distribution by using minimization 
to find the best fit to the data. A two-temperature model i s 
used to describe the star l|Collier-Cameron fc Unruhlll994l ) 
and each element of the star is assigned a spot-filling factor 
/. DoTS finds an optimal solution by minimising while 
maximising the entropy of the image. It is not sufficient to 
simply minimize as a large number of images may fit the 
data, particularly when the errors are large. The maximum 
entropy approach allows a unique image to be returned that 



does not contain information that was not required by the 
data (Skilling & Bryan 1984). DoTS also requires a look-up 
table, containing the shape of the unbroadened spectral line 
at a number of limb angles and at the spot and photospheric 
temperatures. 



4.2 Stellar Parameters 

Many of the parameters of TWA 6 were known from previ- 
ous work and are listed in table O Also listed are our new 
or updated parameters as determined using minimiza- 
tion. Testing the program with artificial data has shown 
that nearly all of the parameters can be derived indepen- 
dently by finding the value which gives the lowest (even 
if the other parameters have not yet been fixed). However, 
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Table 2. As tablc[T]for observations made on the nights of 17/02/06 to 20/02/06. 



U D J ect 


Date 


U i otart 


Exp. 


INO. OI 




Comments 


Conditions 








timc[s] 


Frames 








HR 3037 


17 


17:07:40 


20 


1 


75 


Telluric 


Late start 


TWA 6 


17 


17:10:02 


900 


6 


59 


Target 




GJ 529 


18 


16:53:38 


200 


1 


38 


K4.5 


Storms 


TWA 6 


18 


17:14:09 


900 


4 


67 


Target 




BD -20 4645 


18 


18:23:37 


650 


1 


114 


K6 




HR 3037 


19 


09:55:30 


20 


1 


121 


Telluric 


Cloud/seeing > 2" 


BD -02 801 


19 


10:02:10 


800 


1 


115 


K6 


Cloud/seeing > 2" 


TWA 6 


19 


10:18:29 


900 


3 


42 


Target 1 


3rd aborted 768s 


HR 3037 


19 


15:05:55 


20 


1 


74 


Telluric 


Storm 


TWA 6 


19 


15:08:40 


900 


8 


58 


Target 1 


Some cloud 


TWA 17 


19 


17:19:03 


1200 


1 


46 


Target 2 


Some cloud 


TWA 6 


19 


17:41:09 


900 


4 


54 


Target 1 


Some cloud 


HR 3037 


20 


09:40:29 


20 


1 


110 


Telluric 


Cirrus cloud 


TWA 6 


20 


09:45:46 


900 


14 


38 


Target 1 


Cirrus cloud 


HR 3037 


20 


14:36:28 


20 


1 


97 


Telluric 


Cirrus cloud 


TWA 6 


20 


14:39:22 


900 


10 


51 


Target 1 


Cirrus cloud 


TWA 17 


20 


17:22:40 


1200 


2 


31 


Target 2 


Cirrus cloud 



1,0 



0.8 



0.6 



0,4 



0.2 



0.0 



III 



4 


1 

J 




1 
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-50 50 



Figure 2. Greyscale images of the LSD profiles divided by the 
mean profile for ll-13th Feb (left) and 17-20th Feb (right). Spots 
appear as emission (white). The greyscale is given in the bar at 
the bottom. 



for some parameters the minimization must be carried out 
for two parameters simultaneously, e.g. period and differen- 
tial rotation must be determined together (see section l43] for 
further discussion of differential rotation). Of the previously 
measured values only one changed by a notable amount: 
the projected rotational velocity, v sin i was found to have 
a value of 72 km s~^, somewh at higher than the previously 
measured value of 55 km s~^ (|Webb et al.ll 19991 ). 

Constraining the inclination is slightly more difficult. 



The minimum of the x plot is not as clear as with some of 
the other parameters and the curve is very shallow around 
the minimum x^, found at an inclination of 47°. This dif- 
ficulty in obtai ning the inclinat i on an gle has been noted in 
other work e.g. iKuerster et al.l (|l994l ) where ±10° is given 
as a reasonable range. The errors given in table [3] were de- 
termined by eye from the width of the 'trough' in the vs 
sini plot. 

Using the inclination angle, period and projected equa- 
torial velocity the radius was calculated. The resulting value 
of l.O5J?0 is within one sigma of the value calculated using 
the distance and V - R co l our o f (0.8 ± QA)Rq, using the 
method described in iGravl (|1992| ). It should be noted that 
the distance quoted in table [3] is determined using Hippar- 
cos distances to other members of the TWA and is therefore 
uncertain so the radius determination using the rotational 
velocity and period is the more reliable one. 

Assuming an age of 8 Myr, an effective temperature of 
4000K and a luminosity o f 0.25 the b e st-fitt ing mass 



of the s t ar, fro m both the Alexander et al.! ( 1989) and the 
iKurucj l|l99ll ) opacities (|D'Antona fc Mazzitelll .1994, ). is 
0.7Mq. The luminosity was calculated using the radius and 
effective temperature. 



4.3 Surface Images 

Fig.|3]shows the fits to the profiles for both datasets. For the 
first dataset we achieve a of 1.7, in the second dataset it 
is 2.6. A reduced x^ of 1.0 is not achieved for either dataset. 
The errors estimated in LSD using photon noise statistics 
are often underestimated in the c ase of unpolarised spectra 
(|Donati fc Collier Cameronlll997l ). No scaling of errors was 
carried out in this work, so values greater than 1 are 
expected. The normalised residuals of the fits compared to 
the data were studied and there are no obvious systematic 
effects increasing the or introducing artefacts into the 
maps. The distortions in the profiles are very large, suggest- 
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Figure 3. Fits (lines) and deconvolved profiles (symbols) for TWA 6 during Feb 2006. The fits for ll-lSth are shown in the two boxes 
on the left, and those for the 17-20th are in the two on the right. The numbers on the right of the plots are the rotational phase at 
time of observation. The error bars for the data are drawn. The larger error bars in the last 11 profiles in the second box are due to the 
removal of a feature due to Moonlight in the spectra. 



Table 3. Previo usly-known and updated p aram eters for TWA 
6. Ref erenc es: [1]: IZuckerman fc Son3 1120041. f2]:ICohen fc Kuhil 
19791'). f3]:IWebb et al.l lll999l) . f4l lLawson fc Crausd (|2005|) , [51: 



D'Antona fc Ma2zitelhlll994l') 



Parameter 



Previous 
Value 



New Value 



Reference 



Spectral type 


K7 


K7 


[1] 


Temp. [K] 




4000 ± 200 


[2] 


V magnitude 


12 




[3] 


V - R 


1.19 




[2] 


Distance [pc] 


51 ± 5 




[1] 


vsini [kms"-*^] 


55± 10 


72 ±1 


[3], this work 


Period [day] 


0.54 ±0.01 


0.5409 ± 0.00005 


[4], this work 


Inclination [°] 




47+10 


this work 


Radius [R©] 




1 05+°!'^ 


this work 


Mass [Mq] 




0.7 ±0.2 


[5], this work 


Lum. [Lq] 




0.25 ±0.04 


[3], this work 



ing the star has high spot fiUing factors and sizable tem- 
perature differences on the surface. This may explain the 
previous identification of the TWA 6 as a binary star in 
Ijavawardhana et al.l l|2006l ). as the presence of a large dis- 



tortion at the line-centre may cause it to be mistaken for a 
double-lined binary. 

The maps given by these fits are shown in Figs. |3]and 
[5l In all maps white represents a zero spot-filling factor and 
a temperature of 4000K, black represents complete spot- 
filling and a temperature of 3300K. Both maps show a polar 
spot and features extending out of it down to lower latitudes 
and some isolated spots near the equator. The features are 
consistent between the two maps, and both are consistent 
with Fig. |6] which was obtained by fitting both datasets si- 
multaneously. The structure of the polar spot can be seen 
more clearly in Fig. [T] which shows the star in a fiattened 
polar projection. 

High latitude and polar spots are not observed on the 
Sun, but the presence of a polar spot is commo n to other 
young K stars e.g. LO Peg , (|Barnes et al.ll2005bl ) and V410 
Tau, jSchmidt et al.ll2005l ). In rapidly rotating stars polar 
spots have been explained by Coriolis forces causing mag- 
netic flux tubes to rise almo st parallel to the rotational axis 
(jSchuessler fc Solankil Il992l ). Ma gnetohydrodynamical cal- 
culations in ICaligari et aLllI 19941 ) confirmed that for stars 
with high rotational velocities and deep convection zones 
high spot emergence latitudes are expected. 
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Map for data from 11 -13th February 2006 



Map for all data from 11 -20th Februar/ 2006 
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Figure 4. Doppler image produced from the data from ll-13th 
Feb 2006, shown as a Mercator projection. White represents a zero 
spot-fiUing factor and a temperature of 4000K, black represents 
complete spot-filling and a temperature of 3300K. The greyscale 
is shown on the right hand side, and latitude on the left. Tick- 
marks along the bottom edge of the map show the phase coverage 
of the data. The average filling factor / is 9.5%. 



Mop for data from 17-20th February 200B 
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Figure 5. As Fig.|4]for the second set of data (from 17-20th Feb 
2006). The average filling factor is 10.7%. 
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Figure 6. As Fig. |4] for the combination of both data sets. The 
average filling factor is 9.5%. 




Figure 7. A polar projection of the Doppler map for the full 
dataset for ll-20th February 2006. Dotted circles are at latitude 
intervals of 30° , the solid circle is at the equator. Longitude inter- 
vals of 30° are also shown . The numbers around the edge show 
phases from 0-1, corresponding to those given in Fig|3] Tickmarks 
are placed at phases where observations were made. 



iGranzer et al.] (1200011 carried out simulations of flux 
tube emergence in stars with a variety of rotation rates and 
masses. According to these models, for a TTS with a mass of 
0.7 M© and an angular velocity 50 times that of the Sun, the 
predominate flux tube emergence should be within a latitude 
range 60 — 90°. The lower latitude spots are not predicted. 
In the model for rapidly rotating stars {Q, ^ 10 il©) with 
very small radiative cores the Coriolis force is expected to 
overwhelm the buoyancy forces. The flux tubes detach from 
the overshoot layer altogether and the resulting flux rings 
rise close to the pole, giving rise to little spot emergence be- 
low 60°. TWA 6 has spot emergence at all latitudes above 
the equator, and the spot distribution looks more like that 
for a 1.0 Mq star. However given the uncertainty in its mass 
this is not unreasonable. 



4.4 Photometric Lightcurves 

The Doppler images produced in section 14.31 were used 
to produce monochromatic synthetic lightcurves centred at 
B,V and I. The results were double peaked proflles with the 
secondary peak being about one-eighth the size of the pri- 
ma ry. Photometry is often used as a constraint on spot sizes, 
e.g lVogt et ai] (| 19991 ') . however contemporaneous photome- 
try of TWA 6 does not exist. We can qu alitatively compare 
our li ghtcurves with the photometry in iLawson fc Crausd 
(|2005l ). taken in March 2000. In doing this we are not in- 
tending to directly compare the lightcurves, as several years 
elapsed between the two sets of measurements, only to check 
that our spot maps are plausible in hght of previous results. 

As the central wavelength of our LSD proflles is ~ 
5500A comparison should be made with the V band pho- 
tometry. The amplitude of the variability, 0.1 magnitudes is 
around one-flfth of the value measured bv lLawson fc Crausd 
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l|2005l ') (AV = 0.5). At first sight tliis suggests tliat tlie spot 
coverage wlien our measurements were made was rather dif - 
ferent to when the photometry in iLawson fc Crausd (|2005l l 
was carried out. It should however be noted that it is rather 
difficult to reliably recover the amplitude of the variation us- 
ing spectroscopy alone as spots at low latitudes will be given 
less weight in Doppler i maging, and the am plitude will be 
underestimated, see e.g. lUnruh et al.l (|l995l ). 

The lightcurve amplitude of AV = 0.5 magnitudes is 
large compared to other WTTS, where amplitudes of ~ 0.1 
are more typical. A check should therefore be made to make 
sure our spot/photospheric flux differences are realistic. A 
simple calculation suggests that a spot coverage of approxi- 
mately 40% on one hemisphere (i.e. 20% of the total surface) 
can recreate the observed brightness change. By using DoTS 
to produce synthetic images effects like limb darkening can 
be taken into account. Using our spot temperature of 3300K 
a hghtcurve with an amplitude of 0.5 can be produced with 
a gaussian spot at 30° latitude and a mean spot-filling fac- 
tor of 9%. Such a spot has a full width at half maximum 
(FWHM) of 13°. The mean filling-factor in the images in 
Figs. SI [5] and [6] is 10%. This suggests that there was only 
one major spot group in 2000 whereas in 2006 there were rel- 
atively large mid-latitude spots apart from the main group 
at phase 0.0. Hence it is likely that a combination of un- 
derestimation of the variability using Doppler imaging and 
a change in the spot coverage over the intervening years 
causes the discrepancy. There is certainly a case for further 
photometry on this object, as subsequent sections will show. 

4.5 Differential Rotation 

Figs. |4] and [5] show close agreement between the spots in 
the two datasets. Fig. [8] shows the contours of the second 
map overlaid on the greyscale image of the first. The two 
images in this case were constructed using similar phases in 
order to remove differences that may arise due to the differ- 
ent phase coverage of the two datasets. The alignment of the 
spots suggests that there is no strong differential rotation in 
TWA 6. The main differences are in the spots at lower lati- 
tudes, where the spots in the greyscale seem slightly larger. 
However given the difficulties in determining the exact lati- 
tudinal extent and weight of low-latitude features some dif- 
ferences are to be expected. There is no strong evidence of 
the meridional circulation of the type referred to in section 

To quantify the differential rotation more rigorously the 
two images were cross-correlated at each latitude band with 
a range of longitude shifts. If differential rotation of the type 
observed in the Sun was occurring this would be apparent 
as the band of maximum correlation would curve towards 
negative longitude shifts at higher latitudes. A solar-type 
differential rotation can be described by 

Q.(4>) = f7o - ^^sin(j), (1) 

where fio is the angular velocity of the star at the equator, (j> 
is the latitude and Afi is the difference in angular velocities 
between the equator and the poles. Fig. |9] plots the value of 
the cross correlation in a latitude range — 70° . There is no 
apparent curvature in the position of the maximum corre- 
lation. An upper limit on Afi can be calculated by fitting 
gaussians to each latitude band in the cross correlation. At 




100 200 300 

longitude [deg] 



Figure 8. The greyscale image from first dataset overlaid with 
the contours of the second for the northern hemisphere. The two 
images were reconstructed using similar phases. 




-50 50 

brgitude 5hift ['] 



Figure 9. Cross correlation of the first map with the second. The 
soUd white fine shows the maximum correlation at each latitude. 
The black curves are a contour map of the velocity represented 
by each part of the cross-correlation image. Av is the velocity 
resolution of the instrument (6.7 km s~^). 

each 10° latitude band between — 70° the mean position of 
the gaussian peak is always within one pixel of 0° (one pixel 
corresponds to 2°). At higher latitudes it becomes more diffi- 
cult to cross-correlate the images due to the large polar spot 
and lower velocity ranges at those latitudes. At lower lati- 
tudes the width of the gaussian can provide an upper limit 
to the differential rotation. The mean width of the gaussians 
between — 50° (where most of the isolated spots are) gives 
an upper limit on ASl of 0.04 rad day~^. 

Carrying out a similar cross-correlation, but this time 
comparing longitude bands between the two images, allows 
us to investigate whether the star exhibits differential ro- 
tation along the meridional direction. As before, a gaussian 
was fitted to each band in the correlation image (not shown), 
the mean FWHM of these gaussians gives an upper limit on 
this form of differential rotation of 0.02 rad day^^. 

As a further constraint on the differential rotation DoTS 
was used to explore the parameter space to find the values 
which give the minimum . The relevant parameters which 
should be varied can be found in equation[T] The differential 



Doppler Images and Chromo spheric Variability of TWA 6 9 




1 

11. 612 


11.614 11.616 


1 , , , 
1 1.618 




Q [rad day"'] 







Figure 10. A plot of the value for a variety of angular speeds 
and differential rotation values, where f2o and Af2 are as in equa- 
tion[T] The values are given in the bar at the bottom. The min- 
imum x^, at the '+' in the plot, is at H = 11.616 rad day~^ (cor- 
responding to a period of 0.54091 d) and AQ = 0.0 rad day~^. 
The smallest contour is at a of 2.376, and approximately cor- 
responds to the 3(T confidence level. 



rotation of the star can be parameterised using AfJ, but 
fio must be varied simultaneously, equivalent to varying the 
period and differential rotation. The result of this is shown 
in Fig. 1101 which is a two-dimensional plot of against 
differential rotation and the angular velocity. We find that 
the minimum is at period = 0.54091 day and differential 
rotation = 0.00 rad day~^, further evidence that the star 
does not have differential rotation. 

By drawing contours of Ax^ confidence levels can be 
placed on the values of f2o and An obtained. When us- 
ing x^ minimization to find a parameter value there is a 
67% probability of th e true value lying in the range where 
Ax^ ^ 1 (Press 2 0021 ). The la confidence level therefore lies 
at Ax^ ~ 1. As our minimum is more than 1.0 all the 
values must be divided by the minimum . The surface 
is then multiplied by the number of degrees of freedom in 
the fit, i.e. the number of spectral data points used to create 
the Doppler image. The la confidence level indicates that 
differential rotation lies within the range ±0.003 rad day~^. 



5 BALMER LINE ANALYSIS 

Balmer lines are a useful probe of the circumstellar environ- 
ment of TTS. Ha emission is formed in the stellar chromo- 
sphere, disc and wind when hydrogen atoms are ionised by 
photons from the stellar surface. When a prominence lies in 
front of the stellar disc from our point of view most of the 
Balmer line photons will scattered out of our line-of-sight, 
so prominences in front of the stellar disc appear in absorp- 
tion. When the prominence is off the stellar disc the photons 
will be scattered into our line-of-sight and the prominence 
appears in emission. 

Fig. [11] shows the mean Ha and H/3 profile shapes for 
the spectra in the first and second datasets. The raw profiles 




veloc1ty[km/s] 

Figure 11. Mean Ho profiles (top) and H/3 (bottom). In each 
plot the upper set of lines shows the normalised mean profile for 
the first dataset (solid line) and second dataset (dashed line). The 
lower set of lines are the mean profile with a rotationally broad- 
ened absorption profile subtracted off. The template spectrum 
was rotationally broadened to the vsini of TWA 6 (72 km s"'^), 
shifted to an appropriate radial velocity, scaled and subtracted. 
The telluric lines have been divided out. For Ho the unsubtracted 
and subtracted lines have an equivalent width of SA and 4.2A re- 
spectively, while for H/3 it is 1.2A and l.oA. 



are shown on top and underneath a broadened absorption 
spectrum has been subtracted off. The absorption spectrum 
was created by broadening the normalised template K7 spec- 
trum to a wsin j of 72 km s~^. The telluric lines have been 
removed. 

Both Ha and H/3 profiles have a narrow component with 
a FWHM of 110 km s"^and 113 km s"^ respectively. 
This in good agreement with the expected FWHM of a ro- 
tationally broadened spectral line for a star with a v sin i of 
72 km s"\ which is 112 k m s~^ dC ray 1992) ff a limb dark- 
ening of 0.88 is assumed ||A1-Naimivlll978f ). The subtracted 
Hq profiles have additional broad wings. These are more 
prominent on the blue side. The subtracted H/3 profiles by 
contrast appear to only have a narrow component. 

The presence of a broad component on the Ha emission 
is supported by the calculati on of the variance pro file, shown 
in Fig. [12] As described in jjohns fc Basril (|l995l ) this is is 
calculated using 



(n-1) 



1/2 



(2) 
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Figure 12. The normalised variance profiles for H«. The solid 
line is the variance profiles for the first dataset and the dashed line 
is for the second dataset. This shows variance out to ±250 km s~^ . 
The small peak at ~ —350 km s~^ is due to the presence of an 
absorption line. 



Fig. [12] shows that the Ha emission varies out to 
±250 km s~^, well beyond the wsini of 72 km s-^ A 
similar calculation was carried out for H/3 but this showed 
negligible variability off the stellar disc and has not been 
shown here. 

Figs. [13] and [14] are greyscale images of Ha divided by 
a mean profile, at narrow and wide velocity scales, shown 
with the LSD profiles for comparison. In the left and centre 
images of Fig[l3]there is a strong similarity between the LSD 
and Ha profiles. The bright regions coincide in these images, 
suggesting that the spots on the photosphere correspond 
with active Ha regions. The equivalent H/3 image has similar 
features, although it is noisier and has not been shown. The 
SNR in Ha is lower than for the LSD as the LSD profiles 
have been produced using about 2000 lines; we are thus not 
able to resolve fine detail in the Ha images that we can see 
in the LSD image. 

The correspondence is less clear between phases - 0.4 
in Fig. llSl than in phases 0.4 - 1.0. There is a dark feature at 
around 50 km s~^ in the Ha which is making the contrast 
more difficult to see. In Fig. [14] there is also similarity, be- 
tween phases - 0.4 some of the bright regions coincide. At 
phases 0.4 - 0.6 in the Ha there are dark features masking 
any possible correspondence but there is some correspon- 
dence between phases 0.6 - 1.0. 

Fig. [E] shows the LSD image with the contours of 
the Ha on top. Both images have been smoothed. The 
tick marks show the direction of decreasing intensity. The 
straight horizontal edges in the contour plot are where the 
phase coverage was incomplete. The alignment between the 
bright regions is apparent in the left hand image for phase 
0.4 - 1.0. In the right hand image alignment appears to be 
present between phases 0.6 - 1.0 only. Note that we have as- 
sumed Sun-like rotation, so a phase of 0.4 corresponds to a 
longitude of 216°. Hence the increased Ha emission appears 
to associated with the group of spots between 60 — 220° . 

The behaviour implied by Fig [13] is similar to the 
Sun but has not often been observed in other stars. 
While low-lying plages have bee n inferred previou s ly, e.g . 
iDonati fc Collier CameronI l|l997h . [Fernandez et al.l l|2004h . 



we believe that this is the first observation of such a clear 
alignment of phases and velocities, indicating that the plages 
are overlying the spots. This is in contrast to observations of 
earlier type stars such as AB Dor and Speedy Mic where the 
dominant signal is from higher-lying prominences (known as 
'slingshot prominences'). 

Also shown on the right-hand sides Fig [T3] and [2] is Ha 
variability within a velocity range of ±500 km s-\ There 
are features well outside the v sin i of the star, suggesting 
that slingshot prominences can be seen in emission off the 
stellar disc. In particular there is a bright emission feature in 
Fig |13l that is can be fitted with a sinusoid, shown in Fig. 1161 
In this figure the image has been smoothed and the contrast 
has been enhanced, and a sinusoid has been overplotted. 
The maximum velocity of this sinusoid is 285 km s~^. If we 
assume the prominence is corotating with the star then it is 
lying at a heig ht of about 4fl« . The star has a co rotation ra- 
dius of 2.4flj,. Ijardine fc van Ballegooiier] ()2005l l show that 
a star of this corotation radius can support prominences out 
to 4.8i?, . The model assumes that the prominences are em- 
bedded in the wind rather than an extended corona. Sling- 
shot prominences have not been detected before in T Tauri 
stars - their presence suggests that the star can support 
stable structures beyond corotation. 

Other explanations for this large velocity range should 
also be considered. For example, there may be infalling gas, 
close to the surface but with a maximum velocity component 
along our line of sight of around 285 km s~^. The observed 
velocity will vary sinusoidally as the star rotates. Unless the 
gas is falling directly onto the equator the magnitude of the 
velocity would not be the same at phases separated by 0.5, 
as they appear to be here. However if it was on the equator 
we would not see it for part of the rotation as the inclination 
of the star would cause it to go out of sight. 

The spectra plotted in Fig. [13] were taken over 3 nights, 
corresponding to about 5.5 rotational periods. All the spec- 
tra used in Fig. [13] were taken on either the first or third day 
of observations. Most taken on the third day lie in the phase 
range 0.5 - 0.7 where the brightening is less obvious, which 
implies that the prominence was only present on the first 
night. However there are two spectra from the third night in 
the phase range 0.7 - 1.0 which do support the presence of 
a prominence, in particular the line at phase ~ 0.75 shows 
a brightening at velocities ~ 300 km s~^. The prominence 
may be less apparent at phase 0.5 - 0.7 because it is on the 
far side of the star at that time and is obscured by circum- 
stellar material. 

In Fig. [14] there are features outside v sin i that may 
be evidence of more slingshot prominences, particularly be- 
tween phases 0.6 - 1.0. However they do not seem to be 
present at other phases and there is not enough information 
to justify identifying them as prominences. Also, as they are 
not at the same phases as the prominence in Fig. [T3] it is 
clear that this prominence has not survived the 5 days be- 
tween the two sets of observations, suggesting lifetimes of 
the order of a few days. 

6 CONCLUSIONS 

Spectra of TWA 6 were taken at AAT/UCLES during two 
epochs separated by approximately nine rotational periods. 
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Figure 13. Greyscale images for the first set of data (covering 
11 — 13th Feb 2006). From left to right, the normahsed decon- 
volved profiles from Fig [21 the normalised Ha at velocity range 
±90km s~^ (bin width 4.2 km s^^) and the normalised Ha at 
velocity range itSOOkm (bin widths 17 km s~^). The scale of 
the centre and right images is given by the bar at the bottom, 
while the scale of the left image is as in Fig. [2] The similarity 
between the left and central images should be noted, the bright 
regions correspond, particularly from phase 0.4 - 1.0. In the right- 
hand image there is a sinusoidal feature with an amplitude of 
285 km s~^, see also Fig. 1161 
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Figure 15. LSD profiles shown in greyscale, overlaid with the 
contours of the H« profiles, for ll-13th Feb (left) and 17-20th 
Feb 2006 (right). The LSD and Ha profiles are as in Fig.[T3land 
1 141 but they have been smoothed. 




0.90 0,92 0,94 0,96 0.96 



i C4 1.06 l.OS 1,10 1,12 



:ity [krr 



Figure 14. As Fig. 1131 for the second set of data, covering 17- 
20th Feb 2006. The scale bar refers to the centre and right images, 
the left image is as in Fig. |2] 



Figure 16. Ha variability overplotted with a sinusoid indicating 
the position of a slingshot prominence. The maximum velocity is 
285 km s^^. The vertical lines are at iiisini, i.e. ± 72 km s~^. 



Using these spectra and Doppler imaging two maps of its 
spot distribution were produced. The spot distribution of 
TWA 6 is qualitatively sim ilar to other T Tau ri stars, e.g LO 
Peg, (jBarnes et all l2005bl ) and V410 Tan, (jSchmidt et all 
I2OO5I ). with the predominate spot emergence at high lati- 
tudes and spots all the way to the equator. 

In the Sun photospheric spots are often coincident with 
'plages', active regions in the chromosphere, usually ob- 
served using a Ha filter. The Hq emission of TWA 6 in- 



dicates that the active chromospheric regions are cospatial 
with the photospheric spots - similar to the behaviour of 
the Sun. This supports previous findings that Ha emis- 
sion in young stars incre ases whe n spot coverage i s larger 
iDonati fc CoUier Caniieron (1997, ) , iFernandez et al.l ([2004) ) 
but here we have, for the first time, a strong indication that 
the plages are overlying the spots. 

TWA 6 appears to support slingshot prominences, the 
first such observation in a T Tauri star. At least one promi- 
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nence is observed, lying at a velocity of 285 km s^^ which 
corresponds to a radial distance of 4i?, from the axis. There 
is evidence that the prominence survives for three days after 
the initial detection - suggesting that the structure is stable 
over at least 5.5 rotational periods. The presence of a stable 
slingshot prominence at a height of 3i?* from the surface 
suggests that the star has a large-scale field. Further obser- 
vations of TWA 6 are justified in order to fully characterise 
its prominence system, as well as of other stars of similar 
age and spectral type to determine whether the behaviour 
set out in this paper is typical. 

The star has no differential rotation (upper limit 
0.003 rad day"^). in line with the previously observed trend 
that differential rotation decreases with decreasing effective 
temperature and increasing rotational velocity. Earlier dif- 
ferential rotation measurements have been for older objects, 
here we have extended the observations to younger stars. 
The null-measurement is consistent with expectations for 
the spectral type and angular speed of TWA 6 and indi- 
cates that T Tauri stars exhibit solid-body rotation. 
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